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Abstract 

Background: Wireless energy transfer system is an effective way to solve the visual 
prosthesis energy supply problems, theoretical modeling of the system is the 
prerequisite to do optimal energy transfer system design. 

Methods: On the basis of the ideal model of the wireless energy transfer system, 
according to visual prosthesis application condition, the system modeling is optimized. 
During the optimal modeling, taking planar spiral coils as the coupling devices 
between energy transmitter and receiver, the effect of the parasitic capacitance of the 
transfer coil is considered, and especially the concept of biological capacitance is 
proposed to consider the influence of biological tissue on the energy transfer 
efficiency, resulting in the optimal modeling's more accuracy for the actual application. 

Results: The simulation data of the optimal model in this paper is compared with that 
of the previous ideal model, the results show that under high frequency condition, the 
parasitic capacitance of inductance and biological capacitance considered in the 
optimal model could have great impact on the wireless energy transfer system. The 
further comparison with the experimental data verifies the validity and accuracy of the 
optimal model proposed in this paper. 

Conclusions: The optimal model proposed in this paper has a higher theoretical 
guiding significance for the wireless energy transfer system's further research, and 
provide a more precise model reference for solving the power supply problem in visual 
prosthesis clinical application. 

Keywords: Biological capacitance, Optimal modeling, Visual prosthesis, Wireless 
energy transfer, Planar spiral coils 
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Background 

Visual prosthesis is one of the hot research issues in biomedical field and how to solve 
the energy supply problem of implantable visual prosthesis is particularly important. To 
reduce the patient s pain of multiple surgical implantable battery replacement, in recent 
years wireless energy transfer become scholars accepted method [1-9]. Many scholars 
have done some research on wireless energy transfer of the coupling coils, but the mod- 
els built in these papers are too ideal, for example the internal resistance of power source 
and the parasitic capacitance of the coil are neglected, and most of them are based on air 
medium [1-3]. Part of the literatures have considered the biological tissue influence on 
energy transfer, but they all are based on the analysis of the simulation software such as 
Ansoft-HFSS without establishing more accurate mathematical model on the basis of the 
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physical properties of biological tissue [4-8] . In this paper, we focus on the optimal mod- 
eling of visual prosthesis wireless energy transfer system, During the modeling, on one 
hand, the factors that may affect the wireless energy transfer are fully considered without 
taking many ideal assumptions. On the other hand, the influence of biological tissue, the 
coupling medium between primary coil and secondary coil, is fully taken account. 

Methods 

The ideal modeling of wireless power transfer system 

Figure 1 shows the model of ideal wireless energy transfer (wireless power transfer, WPT) 
system [1], in which the internal resistance of power source r s and the coil parasitic capac- 
itance are ignored. Energy transfer is achieved by the coupling of the primary inductance 
coil L\ and the secondary inductance coil L 2 . Usually the mutual inductance (factor) M 
or the coupling coefficient K\ 2 characterizes the coupling degree of two coils, and the 
formula is K\ 2 = r¥-j- < m Figure 1, V s means the energy source. Z\ means the equiv- 



alent impedance of the primary coil circuit, which is composed by the reactance jcoL\ 
of the primary coil L\ (where co is the energy transfer angular frequency) and equiva- 
lent series resistance r\ , that is Z\ = r\ + j(oL\. Z 2 means the equivalent impedance 
of the secondary coil circuit, which is composed by the reactance jcoL 2 of the secondary 
coil L 2 , the equivalent series resistance r 2 and the load Zi (Zi = Ri + jcoXi), that is 
Z 2 = r 2 + jooL 2 + Z L . 

According to the energy conservation law, the total energy supplied by the energy source 
V s is consumed by r\,r 2 and Ri jointly, and the reactance parts of L\,L 2 and Xl do not 
involved in energy consumption. Therefore, the energy transfer efficiency r] of the wireless 
energy transfer system can be defined as the ratio of the energy consumed by the load and 
the total energy provided by the energy source. That is: 



Figure 1 The model of ideal wireless energy transfer system. Figure 1 shows the model of ideal wireless 
energy transfer (wireless power transfer, WPT) system. 




n 



a) 



In Formula (1), 
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P Rl = t^Ri, which is the energy obtained (consumed) by the load, z 2 is the secondary 
coil loop current. P n = i\r2, which is the energy consumption of the equivalent series 
resistance r 2 in the secondary coil P n = i\r\, which is the energy consumption of the 
equivalent series resistance r\ in the primary coil, i\ is the primary coil loop current. 

h = Z2VS 2 (2) 
ZiZ 2 + (coM) 2 

= -jcoMVs 
2 Z X Z 2 + (coM) 2 

So, the ideal mathematical model of the wireless energy transfer system can be 
summarized as follows: 

coM 2 ^ 

rj — tl (4) 

^| + a>AP% + [(coL 2 + X L ) 2 + (r 2 + R L ) 2 } Jfc 

Modeling optimization of visual prosthesis wireless energy transfer system 

During the modeling optimization in this paper, the impacts of the physical parameters 
in the actual circuits to the wireless energy transfer are fully considered without tak- 
ing ideal assumptions, and especially the impact of biological tissue is taken account. 
Figure 2 shows the equivalent circuit of the perfect visual prosthesis wireless energy trans- 
fer system proposed in this paper. Compared with the ideal model shown in Figure 1, the 
internal resistance of power source r s and the parasitic capacitance of the primary coil C p 
are added in the primary coil circuit impedance Z\, Its worth to note that in this paper 
we will use the vector network analyzer to measure the wireless energy transfer efficiency 
of the system, for the characteristics of the equipment, here we have r s = Rl = 50Q, 
and for the actual application, the load resistance Rl should be the equivalent load of the 
circuit in vivo. In the secondary coil circuit impedance Z 2 , the biological capacitance C p 
proposed in this paper is added. When the secondary coil is implanted into the body, the 
air medium is replaced by the biological tissue medium (the relative dielectric constant of 




Figure 2 The visual prosthesis wireless energy transmission system. Figure 2 shows the equivalent 
circuit of the visual prosthesis wireless energy transfer system proposed in this paper. 
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human biological medium s m >> 1), in this paper the capacitance produced by the bio- 
logical medium is defined as the biological capacitance. In addition, in order to improve 
the energy transfer efficiency, the resonant match capacitances C\ and C2 are added in 
the primary side and the secondary side respectively. 

The modeling of the planar spiral inductance coil L 

The coupling coils used in this system are the planar spiral inductance coils, the shape is 
shown in Figure 3. 

Such Inductance coils are widely used in wireless charger, card reader, IC, ID card, 
which have the following physical properties: 

1. They can be manufactured conveniently and they are economical and durable; 

2. Good flexibility, easy to be implanted. 

Based on large number of experimental data statistics, the planar spiral coil inductance 
empirical formula [10] is adopted here: 



(2r + 2.8d) * 10 5 

From above formula (5), we can get that the parameters affecting the coil inductance 
can be summarized as: 
Coil winding turns n; 

Average coil winding radius r (the average of the outer radius and inner radius, unit: m); 
Coil winding depth d (outer radius minus inner radius, unit: m). 

In this paper, the size of the primary and secondary spiral coil (Figure 3) is the same. 
Number of turns n is 8, average winding radius r is 17.08mm, and winding depth d is 
4.34mm. According to the proposed model we can obtain the planar spiral coil inductance 
that is Li =L 2 = 4.031 /xH. 




Figure 3 The planar spiral inductance coil. The coupling coils used in this system are the planar spiral 
inductance coils. 
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The modeling of the mutual inductance M 

When the central axis of the primary coil and the secondary coil are aligned, according to 
Maxwells equation, the mutual inductance value My between any pair of parallel single- 
turn coils with the radius of r/, rj can be expressed as equation (6) [6]: 

2/x 



My- — 
a 



a 2 



K(a) - E(a) 



(6) 



wherea = 2 /SpT^ ; 
K(a) and E(a) are class I and II of complete elliptic integral respectively; 

f 1 dt 

Class I of elliptic integral expression is: K(a) = I — 

Jo j(l-t 2 )(l-a 2 t 2 ) 



dt 



r 1 Vi - olH 2 

Class II of elliptic integral expression is: E(a) = / — dt 

Jo VI - t 2 

D is the distance between the coils; fi is the transfer medium magnetic permeability 
between the coils. 

In summary, when the turns of the primary coil are n\, the turns of the secondary coils 
are the mutual inductance M between two coils can be expressed as: 

n\ ri2 

M = Y J Y J M ij(n,r j ,D) (7) 

i=l ;'=1 

As the coil is shown in Figure 3, the inner diameter is 14.91mm; the diameter of each 
turn is 0.62mm; the vacuum magnetic permeability is /xn = 47T * 10 -7 H/m; the relative 
magnetic permeability of air and biological tissue is \i r ~ 1. Substituting into this model, 
when the distance between the coils D=lcm, the mutual inductance M = 1.0209/xH; 
when the distance between the coils D=2cm, the mutual inductance M = 0.42134/xH. 

The modeling of the inductance coil high-frequency equivalent series resistance R 

Under high frequencies condition, because the inductance coil resistance is influenced by 
the skin effect, the proximity effect, the eddy current effect and many other effects, the 
impedance is much larger than the intrinsic resistance of coil at DC, Typically, the coil 
loss resistance under high frequency can break down into Ohm loss resistance R 0 and 
radiation loss resistance R r [3,11]. In this paper, the energy transfer frequencies are in 
l-50MHz, resulting in R r << R 0 , so the radiation loss impedance R r can be negligible. 



In above formula, /j,q is the vacuum magnetic permeability; o is the conductivity; / is the 
wire length; r' is the radius of the wire; n is the number of the coil turns; r is the average 
radius of the coils; co is the angular frequency. 

In our system, the coil conductivity a = 5.9 * 10 7 S/m; the average radius r = 17.08mm; 
turns n — 8; the equivalent series resistance is changed with frequency, for example, when 
co = 35.35534rad/s that is/ = 5.63MHz, the coil equivalent series resistance R 0 = 0.2705£1 

The Modeling of the inductance coil parasitic capacitance C p 

Figure 4 is the cross-sectional view between two adjacent turns of the coils. In the figure, 
D c means the bare wire diameter, D 0 means the wire diameter with the insulating layer, 
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Figure 4 The cross-sectional view between two adjacent turns of the coil. 

and x(0) means the air gap between two adjacent wires. The relationship between x(0) 
and other parameters could be given as 



x{6) = D 0 (l-cosO) 



(9) 



The parasitic capacitance C# between two adjacent turns of the coil can be equiv- 
alent to the insulation layer dielectric equivalent capacitance C ttc in series with the 
inter-insulation layer air medium equivalent capacitance C ttg [12]. 

The insulation layer dielectric equivalent capacitance per unit angle can be expressed 
as: 

£ r £ ok Jn 



dCftr — 



(10) 



8 r is the relative dielectric constant of the insulation layer; 

l t is the corresponding effective length of two adjacent turns of the coil; 

The inter-insulation layer air medium equivalent capacitance per unit angle can be 
expressed as: 

UD 0 h 



dC l 



ttg 



£ oz 



£ o: 



-dO 



2x(0) 2(1- cos 0) 
The equivalent capacitance per unit angle can be expressed as: 



dC eq (0) = 



dCttcdCttg _ £oh 1 



dCttc H~ dCt 



ttg 



2 1 + J-ln^ -cos(9 



de 



(ii) 



(12) 



The equivalent capacitance between two adjacent turns of the coil can be expressed as: 



Ct 



= 2 f 

Jo 



2 1 



e r m D c 



-d6 du 



cose 



2s r arctan 



(-l+V3)(2 £r +ln^) 



(l+V3) /ln^(2 £r+ ln%) 



ln^(2, r + ln^) 



(13) 
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By measuring the selected planar spiral coil in this paper we can get the bare wire diam- 
eter D c = 0.575mm, the wire diameter including the insulation layer D 0 = 0.620mm, 
the vacuum dielectric constant So = 8.85 * 10 -12 F/m, the enameled wire insulation 
layer relative dielectric constant s r ^3.5 [13], the innermost effective corresponding 
length between two adjacent turns of the coil l t = 97.3871mm. According to the model, 
the parasitic capacitance between the innermost two turns of the coil can be obtained 
C tt = 9.9399pF; Similarly, the parasitic capacitance between each two adjacent turns of 
the coil can be evaluated. 

The equivalent parasitic capacitance between each two adjacent turns of the multi-turn 
coil also can be expressed as equation (13), and the capacitances between the turns of 
the coil are connected in series, in other words, the lumped equivalent capacitance of the 
n-turn coil equals to the n-1 inter- turn capacitances in series. 

C P = , — i 1 : i— (14) 



C\2 C 2 3 ' ' ' C(„_i)„ 



The Modeling of the Biological Capacitance C p 

The biological capacitance C' p can be equivalent to the insulation layers dielectric equiv- 
alent capacitance C ttc in series with the inter-insulation layer biological tissue dielectric 
equivalent capacitance C ttm (the relative dielectric constant of the biological tissues is 
represented by e m ). The inter-insulation layer biological tissue dielectric equivalent 
capacitance per unit angle could be expressed as: 

kDo h 

dCttm = g 0 g m o = £0£m— —dO (15) 

2x(0) 2(1- cos 0) 

Connect it in series with the insulation layer dielectric equivalent capacitance per unit 
angle C ttc , then the biological capacitance per unit angle can be obtained: 

. , dCff r dCff m £()lf 1 , , 

dC' (&) = -, , = — ; =; ; dO 16 

e l y dC ttc + dC ttm 2 J, + J_i n ^_ J_ cos(? 
The biological capacitance between two adjacent turns of the coil is expressed as: 

(-1+V3) £w lngj 



Jo 



dC' eq (6) = e 0 l t - 



(2s r + s m ln^)e ra ln^ 



(17) 



The value of relative dielectric constant of the biological tissues is different for each one 
and it could be decrease slowly with the increase of the frequency when the frequency 
arrives several MHz [13,14]. In this paper we take the relative dielectric constant s m ~54 
[14], the biological capacitance of the innermost two turns C^=20.192pF. Similarly, 
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So, the modeling built in this paper can be available not only for the general air medium 
wireless energy transfer system but also for the biological tissue medium wireless energy 
transfer system. Formula (14) is used in general air medium and formula (18) is used in 
biological tissue medium. 

Results 

Simulation and analysis 

Using the optimal visual prosthesis energy transfer model proposed in this paper, we get 
the energy transfer efficiency at different frequency point (the distance of the primary coil 
and the secondary coil is 1cm), the "Frequency - Efficiency" curves about the ideal model, 
the optimal model of air medium and biological medium are shown in Figure 5. 

The three curves in Figure 5 show energy transfer efficiency versus the frequency by 
using the resonance method. Among them, the first one is the curve of the ideal model, 
the second one is the curve of the optimal model with air medium, and the third curve is 
the optimal model with biological tissue medium. 

From Figure 5 we can find that at the low frequencies, the three curves are almost 
overlapped, while with the increase of frequency, the difference between each other is 
increasing. The reason is that when the coil inductance is fixed, the resonance frequency 
in the low frequency band needs the larger coupling capacitance, the influence of the air 
medium inductance parasitic capacitance and the biological capacitance is negligible, so, 
the energy transfer efficiency of the model before and after optimization is basically equal. 
However, with the resonant frequency increases, the required matching capacitance grad- 
ually reduces, when it reduces to a certain extent, the influence of the inductance parasitic 
capacitance and its biological capacitance becomes significant, so the trend of three 
curves become discrete. 
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Figure 5 The comparison chart of the energy transfer efficiency simulation (the distance between the 
coils D = 1cm). The "Frequency - Efficiency" curves about the ideal model, the optimal model of air medium 
and biological medium are shown in Figure 5. 
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Experimental verification 

The photos of the experimental devices for measurement is shown in Figure 6, the energy 
transfer efficiency is measured by the vector network analyzer Agilent E5071C. The influ- 
ence of the human biological tissue on the energy transfer efficiency is simulated by using 
the fresh lean from the intraday slaughtered pig wrapping tightly around the secondary 
coil 

Figure 7 is the situation that the energy transfer efficiency varies with the frequency. In 
order to reflect the versatility of this model, in Figure 7 the distance between the primary 
coil and the secondary coil is changed into 2cm. Figure 7 A is the comparison chart of the 
ideal model, the optimal model in the air medium and the measured data. It can be seen 
from Figure 7 A that the modeling simulation results (curve 2) in the air medium is more 
consistent with the experimental data (the point curve 3). The ideal model (curve 1) which 
ignores the parasitic capacitance is basically consistent with the measured data at the 
low frequencies, but the extent of error in the high frequency band cannot be tolerated, 
therefore the optimization modeling theory of this paper is closer to the actual situation. 
Figure 7B is the comparison chart of the simulation data in the biological media and the 
measured data, we can see that the trend of the two curves basically coincide, the extent 
of error range is within the reasonable range, which verifies the correctness of the model. 

The inductance values of the planar spiral inductance coils used in this paper are 
L\ = Z,2 = 4.031/xH, when the match capacitances are 22pF and 47pF, the resonant 
frequencies are 16.97MHz and 11.61MHz respectively. Figure 8 shows that the energy 
transfer efficiency in the air medium and in the biological tissue medium varies with the 
transfer distances at two different frequencies 16.97MHz and 11.61MHz. From Figure 8 
we can find that the air medium modeling highly coincides with the actual measurement 
data, the biological tissue medium modeling also basically coincides with the measured 
data within the distance of l-2cm. According to clinical medicine experience, the dis- 
tance between the secondary visual prosthesis reception coil which implant into the 
patients body and the primary coil outside the body is usually l-2cm. Therefore, this 
model meets the requirements of the implantation depth range for the actual clinical 
application. 

Conclusions 

Based on the physical characteristics of the biological tissue, optimization modeling 
of the visual prosthesis wireless energy transfer system is established in this paper. 




Figure 6 The photos of the experimental device for measurement. A-The experimental instrument and 
the vector network analyzer Agilent E5071 C; B-The picture of the secondary coil embedded by biological 
tissue. 
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Figure 7 The comparison chart of the simulation and the experimental measurement (D=2cm). A-The 

comparison of measurement data in air medium with simulation data using ideal model and optimal model; 
B-The comparison of measurement data in biological medium with simulation data using optimal model 
simulation. 



During the optimization modeling, the influence of the parasitic capacitance of the 
coil and the biological capacitance to the energy transfer efficiency is considered fully, 
which greatly improves the matching degree of the theoretical modeling and the mea- 
sured data. The optimal modeling has a higher guiding significance for the actual 
system design. 





D 



Distance(mm) 



Figure 8 The diagram of the energy transmission efficiency vs. the distance between the coils. A-The 

curves in the air medium at 1 6.97MHz; B-The curves in the biological medium at 1 6.97MHz; C-The curves in 
the air medium at 1 1 .61 MHz; D-The curves in the biological medium at 1 1 .61 MHz. 
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